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DEVELOPMENT OP MICHELL MINIMUM WEIGHT STRUCTURES 


By Dhanjoo N. Ghista and Meyer M. Resnikoff 
Ames Research Center 


SUMMARY 


The Michell strain criterion is used to develop minimum weight structures 
to equilibrate some interesting force systems of significance to aeronautics 
and the space technology. Michell prescribed a strain criterion for the 
members of a minimum volume structure (i.e., a structure having the minimum 
amount of material) to equilibrate a given force system. 

The research work contributory to this paper is divided into three inde- 
pendent sections. In the first section, a theorem is presented on the super- 
position of Michell structures and applied to develop optimum structures for 
two basic force systems. With this theorem a force system that does not lend 
itself to a unique solution (i.e., a unique optimum structure) may be treated 
by resolving it into two or more systems, determining the solutions to the 
resolutes, and then superposing them to give the solution for the total force 
system. The theorem shows that the solution obtained by this superposition 
method satisfies the Michell optimum structure. 

In the second section a solution is obtained to the governing equations 
developed by Hemp for the form of an optimum structure satisfying the Michell 
strain criterion condition in two dimensions. This solution gives a unique 
system of coordinate curves for the member layout of the Michell structure. 
This system can be employed to offer optimum structures for a variety of 
force systems acting on plane rectangular and curvilinear domains . 

In the third section two well-known orthogonal systems of curves - the 
dipolar and the spherical - are selected for member layout of Michell struc- 
tures. The volume of a Michell structure is a function of the displacements 
of the points of application of the applied forces; hence for the dipolar and 
spherical coordinate systems, the generalized displacements are obtained for 
the Michell strain pattern. Some interesting loading systems result when 
the two coordinate systems are used in developing Michell optimum structures. 
The dipolar system gives rise to a structure which could conceivably be 
employed for an entry body that had a displaced center of gravity with respect 
to the geometric center. Another application of this system is a structure 
which distributes a concentrated load over a wide bearing area, thereby reduc- 
ing the possibility of large deformations of the material. An application of 
the spherical system gives rise to a spherical prestressed structure, under a 
suitable force system. 



IOTRODUCTION 


In aeronautics structures are needed that can withstand a given load 
system and that have a minimum volume, of material. For a given load system 
this objective is more easily achieved if the probability of failure is 
allowed to be small instead of zero (ref. 1 ) . However , the development and 
application of the theory of optimum structures, which is the subject of this 
paper, is strictly a problem in structural mechanics; hence, in this case, 
the idealized version of failure (of probability zero) is adopted. 

For research in the field of structural optimization, the structural 
design is usually represented by design parameters of the appropriate shape 
and size, the merit function is expressed ip terms of these parameters and 
minimized with respect to them, within limits of the constraints, by using 
nonlinear programming techniques (see ref. 2) . On the other hand, a theory 
based on displacement or strain criterion could be used to find the structure 
with the absolute minimum volume for the given load system. The form of the 
structure could be developed as a consequence. Such an approach gives the 
absolute minimum volume structure for the given load system. Micheles paper 
(ref. 3 ) inspired research along these lines. For a given system of self- 
equilibrating forces, Michell has laid down a deformation criterion which if 
satisfied by a framework will give a structure of a minimum volume, that is, 
of a minimum amount of material. Cox (ref. k) drew the attention of engineers 
to Michell T s results and, by contributing to the application of these results, 
championed the cause of optimum structural design. Hemp (ref. 5 ) has incor- 
porated the Michell criterion of strain into a two-dimensional theory of 
optimum frameworks and has developed some special forms. Further, Michell 1 s 
theorem has been applied to develop optimum frameworks for some useful 
practical load systems (see refs. 6, 7, 8, and 9)* 

The present paper is presented in three independent parts. In the first 
part a theorem is presented on the superposition of Michell structures. This 
theorem is then used to develop structures for two force systems. In the 
second part, solutions are presented to the governing equations developed by 
Hemp for the form of the minimum volume framework. The resulting structure 
is then developed, that is, the coordinate curves along which the members lie 
are laid out. It is then shown how certain arbitrary functions of the solu- 
tion can be adjusted to match a given external force system. In the third 
part, two systems of well-known orthogonal coordinate curves, namely, the 
dipolar and the spherical coordinate systems, are selected for the layout of 
members of the Michell optimum structures. The corresponding strain- 
displacement equations incorporating the Michell strain criterion are solved 
to give the displacements in terms of corresponding curvilinear coordinates. 
Some force systems are illustrated as applications of the two coordinate 
systems in the development of Michell optimum structures. 
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SYMBOLS 


A,B 

dsi,ds 2 

e 

e aa^ e P|B > e yy 
e <x$’ e $y> e y<x 

Pi 

Ti,T 2 

■ti.>t 2 

(u,v) 

(u a ,u p , u 7 ) 

V* 

a>P>7 

°c 

a t 


unit arc lengths along a and (B members 
arc lengths along a and (3 members 

magnitude of strain in any member of a Michell optimum 
structure 

normal strains in the curvilinear system (a ,(3,7) 
shear strains in the curvilinear system 
forces acting at the points of application r^_ 
forces per unit length parallel to the a and [3 members 
thickness of the a ,|3 members 

displacements corresponding to the curvilinear system (a, 3) 

displacements corresponding to the curvilinear system (00,3,7) 

volume of a Michell optimum structure 

curvilinear coordinates 

allowable stress in compression 

allowable stress in tension 

angles made by the a and (3 curves with the x axis 


MICHELL CRITERION OF STRAIN 


The Michell criterion of strain for a minimum volume framework to 
equilibrate a given force system states that the members of the structure must 
all be strained by the same amount (e), the sign depending on the sign of the 
axial stress carried by the member. 

The volume V of the minimum volume framework to carry a system of 
forces with allowable stresses a c and in compression and tension is 

given by (see ref. 3) 
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(i) 


V 


* 


g t + g c 
20 t a c 



i=x 



i=i 


where r ± denotes the points of applications of the forces F^ and v ± denote 
the displacements of these points. 


This deformation prescribed by Michell imposes certain restrictions upon 
the layout of the members of a Michell structure. At a node of this frame- 
work, the members follow the principal directions of strain. If members 
carrying loads of the same sign meet at a node, there is no restriction on 
their layout since for this case the Michell deformation Is a pure dilatation 
or contraction and the principal axes of strain are indeterminate. If at a 
node, members carrying loads of opposite sign meet, these members must be at 
right angles . The members of this class of frames form curves of orthogonal 
systems . 


SUPERPOSITION OF MICHELL STRUCTURES 


In this section the indeterminacy of Michell structures is discussed. 

Then the utility of superposing Michell structures is presented. A theorem is 
presented on the superposition of Michell structures which relates the volume 
of the superposed structures for the corresponding force systems to the volume 
of the structure for the resultant of the force systems. The usefulness of 
the superposition theorem is demonstrated by means of two examples. 


Indeterminacy of Michell Structures 

In the development of a Michell structure to equilibrate a given force 
system, the geometry of the layout of the members is also unknown in addition 
to the forces in the members. The problem of determining the sizes of the 
members, for an optimum member layout, to equilibrate a given force system 
involves the solution of the following equations: (i) equilibrium equations 

obtained by considering the equation of an element da by dp, where a,p are 
the curvilinear coordinates; (ii) equations of compatibility of strain in the 
curvilinear coordinate system; and (iii) equilibrium conditions along the 
boundary. Once the members f sizes are known, the forces in the members are 
obtained by specifying the allowable stress in the members . 

It is seen that, for Michell structures, the areas of the members and 
hence the forces in the members for an allowable stress are obtained from the 
deformation pattern or kinematics. For a statically determinate structure, 
the determination of member forces for a given external force system requires 
the solution of only the equilibrium equations. A Michell structure is not 
statically determinate since, in addition to equilibrium equations, the 
deformation equations are needed to determine the forces in the members. 
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Superposition 


The strain pattern prescribed by Michell for a minimum volume framework 
is such that solutions satisfying the Michell criterion of strain exist only 
for a few special cases. It often happens that it is not possible to find an 
optimum framework to equilibrate a force system. However, if the force system 
were resolved into component force systems., optimum frameworks could be deter- 
mined for the component force systems . Suppose a composite structure is 
obtained by the superposition of the constituent Michell structures while con- 
tinuity of displacements between the superposed structures is insured. To 
determine whether this composite structure can be used as the optimum struc- 
ture for the given loading in the absence of a unique solution for that load 
system, a theorem is proposed which states that the volume of the composite 
structure is the same as the volume of the unique structure for the given load 
system. This theorem justifies the use of the composite structure as the 
optimum structure for the load system considered. The theorem will now be 
enunciated, proved, and applied to two force systems. 

Theorem : Let F 1; j_, i = 1, 2, 3, denote a self-equilibrating system of 

forces acting at points r^, i =1, 2, 3, and let Vi denote the volume of 
the Michell structure Sj_ which equilibrates the force system . Let 
F 2 i, j = 1, 2, 4, denote another self -equilibrating force system acting at 
points r j , j = 1, 2, 4, and V 2 denote the volume of the Michell structure 

5 2 which equilibrates the force system F 2 j • _Let F 3 ^^ k = 1, 2, 3, 4, 
represent the resultant of the force systems Fii and F 2 j acting at points 
r k , k = l, 2, 3, 4, and let V 3 denote the volume of the Michell structure 

5 3 which equilibrates the force system F 3 j^. If a structure be obtained by 
the superposition of the two Michell structures Sx and S 2 such that there is 
a compatibility of strain between the two superposed structures (i.e., there 
is a continuity of displacements at nodes common to both the structures), then 
the volume of this composite structure will be Vx + V 2 . It is proposed, 
then, that V 3 = Vx + V 2 . 

Proof : Consider a domain D containing the points r 1; r 2 , r 3 , r 4 . In 

this domain an orthogonal coordinate system of curves (or systems of curves, 
provided there is continuity of displacement along the lines of junction) 
links points r x , r 2 , r 3 , r 4 . The Michell strain pattern is now imposed on 
this coordinate system so that the strains along the orthogonal directions are 
±e and the shear strain is zero. Let the displacements at the points r 1; 
r 2 , r 3 , r 4 be y 1; v 2 , . v 3 , v 4 , respectively. 

how at the points r^, i = 1, 2, 3; the system Fxi is made to act. The 
corresponding Michell structure Sx; to equilibrate the system Fxi; will con- 
sist only of those members of the coordinate system that transmit forces; the 
remaining members of the coordinate system can be said to have zero area and, 
hence, do not form part of the structure Sx* Then the volume Vx of Sx 
is (from Eq. (l) ) , 
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Vi = 


g t 4 g c 
2 g t g c 
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Fli 



g t - g c 
2 g t g c 
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Fii * r i 


( 2 ) 


Next, at points r^, j = 1, 2, ^ , the system F 2 j is made to act. The 
corresponding Michell structure S 2 , to equilibrate the system F 2 j , will 
again consist only of those members that transmit forces. The volume V 2 of 
the structure S 2 is given by 


V 2 


g t + g c 

2 g t g c 


I 

«j=l, 2 , 4 



. (I£) . g t , -„ g c 
V e / 2 g t g c 


I 



0 = 1 , 4 


(3) 


The system F 3 k, k = 1, 2, 3, 4, is now made to act at the points r k . 

The volume V 3 of the corresponding Michell structure S 3 , to equilibrate the 
system F 3 k> is given by 


V 3 = 


g t + g c 
2 g t g c 


I 


'vb 

Fak * i-rJ “ 


k=l,2,3,4 


g t - g c 
2 g t g c 


Fak * ik 




k=l,2,3,4 


Wow since F 3 ]j is the resultant of systems Fii and F 2 j , it follows that 

F 3k = Fii + F 2 j 

that is , ~ 

F 31 = Fn + F 2 i 
F 32 = Fx 2 + F 22 

F 33 = F i 3 

F 3 4 = F 2 4 

Structures Si and S 2 _have common points ri and r 2 at which they have equal 
displacements and v 2 . Hence the structures Si and S 2 can be superposed. 
The resulting composite structure will have the volume (Vi + V 2 ). From equa- 
tions ( 2 ), ( 3 ), (4), and ( 5 ), it is seen that the volume V 3 of the Michell 
structure S 3 which equilibrates the resultant of the force systems Fii an< ^ 
F 2 j is equal to the volume (Vi + V 2 ) of the composite structure obtained by 
superposing the two Michell structures S x and S 2 corresponding to the load 
systems Fii and F 2 j ; that is , 

V 3 = Vi + V 2 (7) 


(5) 


( 6 ) 
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Application of the Superposition Theorem 

The theorem is now used to develop optimum structures for the load 
systems of figures l(a) and 2(a) . The load system of figure l(a) , represent- 
ing a pure shear beam-type loading, is resolved into two constitutive systems, 
as shown in figure l(b) . Each system lends itself to a solution, in the sense 
that there is a definite Michell structure to equilibrate each system. These 
structures are shown independently in figure l(c) and are combined as shown in 
figure l(d) to give the optimum structure for the load system under 
consideration. 


p 



M = Pt /2 



(a) Force system 1 (pure shear beam -type loading) . 


P 


m 



•t/Z 


i *2 

f 


M = PI / 2 


(b) Component systems for above force system. 




M 


(c) Michell structures for component systems. 



(d) Composite structure (optimum structure 
for the above force systems) . 


Figure 1.- Michell structure for a pure shear 
beam-type loading. 
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The volume of each Individual unit is (see ref. 9 ) 


V c 





jtM 

2 




( 8 ) 


P 2P P 



2P 2P 


(a) Force system 2 (simple overhang beam-type 
loading ) . 


P P 2P P P 



2P P P 2P 


(b) Component systems for above force system. 


The volume of the Michell 
structure (fig. 1(d)) to carry the 
load system of figure 1(a) is, 
then, 

V* = 2V C 





(d) Composite structure (optimum structure 
for above force system). 


The load system of fig- 
ure 2(a) is resolved as shown in 
figure 2(b) . The Michell struc- 
ture for each component load 
system is shown in figure 2(c) . 
The units are combined as shown 
in figure 2(d) . The volume of 
each unit is (see ref. 9) 


T = 3 ?) < 10 > 


The volume of the Michell struc- 
ture (fig. 2(d)) to carry the load 
system of figure 2(a) is 


Figure 2.- Michell structure for an overhang 
beam-type loading. 
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v* = 3V C 


= — PZ 
2 




( 11 ) 


THE DEVELOPMENT OF A SYSTEM OF COORDINATE CURVES FOR 
MEMBER LAYOUT OF A MICHELL OPTIMUM STRUCTURE 


The governing equations for the form of the orthogonal coordinate curves 
and the equations of equilibrium of a curvilinear element are presented in 
this section. These equations have been derived by Hemp (ref. 5); however , 
their derivation will be briefly represented to maintain continuity of thought. 
A solution is then obtained to Hemp’s governing equations which dictate the 
form of the orthogonal coordinate system that satisfies the condition of com- 
patibility for the Michell criterion of strain in two dimensions. This solu- 
tion gives the layout of coordinate curves which constitute the lines of 
principal strain along which the members of the corresponding Michell struc- 
ture lie. These curves are then enclosed within a suitable boundary by match- 
ing, at the boundary, the forces in the members of the structure and the 
external forces acting on the structure (one gets the nature of external 
forces on the boundary of the structure) . Indirectly, then, one obtains the 
force system which the developed Michell structure (the members lie along the 
coordinate curves and are represented by the solution obtained from Hemp’s 
governing equations for the form of Michell structures) equilibrates. 

Inversely, one can say that the solution (i.e., a Michell structure) has been 
obtained for this force system. 


Equations Governing the Form of the Coordinate Curves 

The outline of the development of Hemp’s governing equations for the 
form of the Michell optimum structure is now presented. (For details see 
section 3 of ref. 5.) Let x - x(a,3), y - y(a,3) be a set of orthogonal 
coordinate curves representing lines of principal strain. The parametric pair 
(a, 3) represents curvilinear coordinates in the plane of the rectangular 
Cartesian coordinate system. The functions x(a,3) and y(a,3) are continuous 
and have first and second derivatives with respect to a and 3- Along 
a-coordinate curves, a varies and 3 is constant; along 3 -coordinate 
curves, 3 varies and a is constant. Positive directions along both of 
these curves are those along which a and 3 are increasing. The arc lengths 
dsx and ds 2 along a and 3 curves are given by 

ds x - A(a,3)da, ds 2 = B(a,3)d3 (12) 


9 



where 


and 


A = 







(13) 


Let ti and ty 2 represent the angles made by the positive tangents to 
a and p curves, respectively, with the x axis such that 

w = i)r 2 _ \jfa_ (lh) 

It follows that 


(15) 


Now 

cos <o = cos(t 2 - ti) 


COS = 


sin i|fi = 


dx 

dsx 

1 3x 

_ A Sa ’ 

i dx 

°° s = SJ 

dsi 

= A 3a ' 

sin „ 2 = §L_ 
ds 2 


1 3x 
' ’ B 5F 


B 3p 


1 /3x 3x 3y 3y \ 

AB \5a 5)3 Spy 


(16) 


From equations (13) and (l6), the derivatives 3^1 /3a and 3c|/ 2 /3p are given 
by 1 


3^i _ 1 3A 3 — ^ 

da B sm co 3p oa ' ' 


3\|/p _ 1 

3p A sin co 


3B 3 

SZ - 5p (A 005 


(17) 


Also, from equations (lk) and (17) the following relation can be shown to 
exist between A, B,. and TZ: 


a f 1 

da |A sin co 


dB d / A — V 

ST - 5F (A 003 ^ 


3u>[ 3 f 1 

S3 T + 3P 1b- sin a) 


5 A d / _ — v 

Sp - < B cos ^ 


= 0 

(18) 


For co = it/2, this relation becomes 


1 For details refer to appendix A of reference 5 • 
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(19) 


d /I dB\ d /l 

5S Va 5^ + 5p v? apy 


= 0 


In the curvilinear coordinate system, the deformation is characterized by the 
strains e^, ego, e a0* The expressions for the strains are given by (ref. 5, 

p- 59) 


_ 1 Su _v_ dA 

~ A Sa AB 5ff 


i av u aB 

e PP B 5p AB 5x 






e a,p 


1 b a f£\ a a AA 

2 [A 5a \B/ B W \A/J 


(20) 


where 


u and v are the displacements along 


a and 0 curves, respectively. 


In the deformed state, the arc lengths dsi and ds 2 are increased by 
1 + eqq and 1 + e™, respectively, and the angle (H is deformed from xc/2 
to (jr/2) - aa0. Replace A, B, and HI by A(l + e^) , B(l + ej, and 
(xr/2) - 2a0, respectively, in equation (l8); if the resulting equation is 
developed correctly to the first order of strain, the equation of compatibil- 
ity of strain is obtained as follows : 


a /b Se pfR . a fk ^eocR ^ 2 e a|3 a Ti aB / 

3 ^ Va —) + Vb ~w) " alap " 5 a L a (eaa ~ e PP } 


a i_ aA / 
ap b ap ^ Saa 



0 a /i aA y 

2 SSTVa §P e aPy> 


5 A 5b 


2 £ 


50 \B 5a ^0 


= 0 


( 21 ) 


Now the Michell criterion of strain, for members carrying strains of opposite 
sign, is defined by 

e aa = e > e pp = “ e > e ap = ® (22) 

For this strain condition, the equation of compatibility of strain becomes 


a /i aB'N a /i aA 

5a \A Sa/ 5F \B apy 


(23) 


It follows from equations (19) and (23) that 
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Thus, 


1 5A 

b 5p 


-F[(a) 




1 8B 
A 5a 


* FaO) 


( 25 ) 


T t 

where F x (a) , F 2 (p) are arbitrary functions and F x (a), F 2 (p) are their deriv- 
atives. The following transformation is now introduced: 


a - cpi(a) , p = cp 2 (p) (26) 

where cp x , cp 2 , cp x , and cp 2 are continuous functions. Further, cpi and cp 2 are 
chosen so that 


9l (a) - Fl 1 (±50 , cp 2 (p) - F 2 1 (±p) 


(27) 


where F^ 1 , F^ 1 are the inverse functions of F x , F 2 and the upper or lower 
sign is taken accordingly as F X F 2 is positive or negative. The transform 
of equations (25) can now be written in terms of a and p as follows: 


1 5A = 1 5B = 

B 5F ’ A 5a 


(28) 


These are the governing equations for the arc lengths A and B of the system 
defined by equations (12) . 

Now, the angles i[r x and t 2 , defining the directions of a and P curves, 
are obtained from equations (l4) with m = n / 2, (l7) and (28) as follows: 

i x =a+p, = |- + a + P (29) 

In order to select a set of orthogonal coordinate curves (l2) that satisfy 
the strain conditions (22), the corresponding rates of changes of arc lengths 
A and B have to satisfy equations (28) . When A, B have been found, equa- 
tions (12), along with equations (29), can be employed to determine the 
equations of the lines of principal strains. 


Equations of Equilibrium 

The equations derived by Hemp (appendix B, ref. 5) for the equilibrium 
of a curvilinear element da by dp are now presented. Let t x and t 2 
denote the thicknesses of the a and p members . Also, let T x , T 2 be the 



forces per unit length parallel to the a and P curves and S be the shear 
per unit length- From the equilibrium of a curvilinear element da by dp, the 
following equations are obtained: 


^(HTJ + + || s - if T 2 » 0 

ST CBS) + ^(AT a ) -!jK + 5§S-0 


(30) 


For Michell criterion, the shear stress equals zero. Further, from equa- 
tions ( 22 ) it follows that 

Ti = h t t i > t 2 = cr c t 2 (31) 


and hence the equations of equilibrium can be written as follows: 


zk ^ Bcr t t i) - 5I (^2) - 0 

^ ( AtJ C t 2) “ || ( ff t t l) = 0 


(32) 


Determination of a System of Coordinate Curves and the Development 
of the Corresponding Michell Optimum Framework 

A solution is now sought to the governing equations ( 28 ) . Hemp has 
obtained a few solutions in the form of certain functions for the arc lengths 
A and B and has presented some layouts (i-e-, coordinate curves) for the 
members of the Michell optimum structures. Herein, solutions to the govern- 
ing equations are obtained in terms of different functions for A and B, which 
then give us a different layout (i.e., a set of coordinate curves) for the 
members of a Michell structure. This layout, then, gives Michell structures 
which can equilibrate new force systems; thus, solutions (Michell structures) 
are obtained for these new systems. To determine a system of coordinate 
curves that satisfy the Michell criterion of strain (22), equations ( 28 ) have 
to be solved- The following functions for the unit arc lengths A and B 
satisfy equations ( 28 ): 


A = k cos(a + p) 

B =• k sin(a + p) 

Equations (15), (29) > and ( 33 ) then give 

x = J cos dsi =■ 4 J cos 2 (a + p)da 

— sin 2 (a + p) + 2 a + Ci(p) 


(33) 
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where Ci(3) is an arbitrary function of p Also, 

x = J cos i|/ 2 ds 2 = -4J* sin 2 (a + p)dp 

=■ -2P + sin 2(a + p) + C 2 (a) 


where C 2 (a) is an arbitrary function of a. On identifying the two expres- 
sions for x, it is seen that Ci(p) is indeed -2P and C 2 (a) is 2a. Hence, 
the expression for x is given by 

x =■ 2(a - p) + sin 2(a + p) ( 3 ^) 


Similarly, 


y - -cos 2(a + p) 


(35) 


Now, eliminating a and p , the Cartesian forms of the a and 3 curves 
senting lines of principal strain) for -1 < y < 1 are, respectively, 


x =■ -43 + cos" 1 (-y) + s/l - y 2 , 
x — 4c& — cos" 1 (-y) + n/i - y 2 , 


a curves 
3 curves 


(repre- 


(36) 


The above equations represent the coordinate curves for the layout of the 
Michell structure. The curves (illustrated in fig* 3) are enclosed in a 
rectangular domain -1 < y < 1, -rt < x < , and can be used to solve plane 
stress problems for rectangular strips loaded with suitable self -equilibrating 
forces. The a and 3 curves are in tension and compression, respectively. 
These tension and compression members are indicated in figure 3 . From equa- 
tions (15) the angles and a|/ 2 made by the a and p curves with the x 
axis are given by 



Figure 3.- Coordinate curves forming a Michell structure within a rectangular domain. 



I 


COS tl 


-J 


l - y 

2 


COS ^2 


. l l t z 

N / 2 ' 


sin \|/ x = / 
sin ^2 


+ 




(37) 


The areas of the members (i.e., thicknesses t x and t 2 of the a and 3 
members) of the structure are now determined from the equilibrium of a curvi- 
linear element of the structure (eqs. ( 32 )) and from the force equilibrium at 
the boundary. If a-^ - a - -cr c in the equilibrium equations (32) and equa- 
tion ( 34 ) is used, the equilibrium equations can be written: 


<ti + t.) + 2(1 + y) 57 1 * 0 


(ti + t 2 ) - 2 (l - y) 


dtg 

By 


i 


1 + y 

2 


= 0 


( 38 ) 


The solution to the above equations is given by 

ti = f(x) + g(x) 


t 2 » - g(x) 


(39) 


where f(x) and g(x) are arbitrary functions of x. It is seen that y = 1 
presents a singularity. The functions f(x) and g(x) are determined from the 
total force equilibrium on the boundary. 

The variation of the forces in the structure for | I ™ I a t I “ a i s 
given as: (i) The tangential force per unit length (Txy) is proportional to 


t x cos \|/ x + t 2 sin \|/ x 


- w 


1 


(40) 


where t x and t 2 are given by equations ( 39 )* (ii) the normal force per unit 
length parallel to the x axis (T x ) is proportional to 


ti cos - t 2 sin ii =• ti J 1 — - t 2 J 1 g y 


(41) 


(iii) the normal force per unit length parallel to the y axis (Ty) is pro- 
portional to 

tj 1 -* 


ti sin ti - t 2 cos \|r 


X-txT 1 


(42) 

15 


I 




Although the external forces 
carried by the structure should be 
self -equilibrating, their distribu- 
tion is governed by the distribu- 
tion 

( 4 l), 

The above development shows 
how the system of coordinate curves 
under consideration develops the 
Michell optimum structure for a 
rectangular domain under a system 
of external forces. These coordi- 
nate curves could also be enclosed 
between an a curve (say (3 - 0 ) , 
a 3 curve (say a = it/ 4 ); and a 
y ordinate (say x - 0 ) . When 
thus bounded, as shown in figure 4 ; 
a concentrated force P at 
3-0 (i.e.; at x - 1 + (it/ 2), 
y =■ 0) can be equilibrated by means 
of normal and tangential forces 
distributed along the y ordinate 
x - 0 . 

Alternatively; the following 
solutions to equations (29) can be 
had: 

A - 4 sin(a + 3) 

B - -4 cos(a + 3) 

The equations for these curves in 
a rectangular Cartesian coordinate 
system are given by 

y - 4a - cos" 1 ( -x) 
y - 43 + cos“ 1 (-x) 


The corresponding curves are illus- 
trated in figure 5 - In nature; 
this system of curves is the same 
as the one considered earlier 
(given by eqs. (33)) except for a 
90° clockwise rotation about the 
origin . 



(43a) 



0 (42)T Ty ' Txy {eqS ' ' 


Figure 5.- Coordinate curves forming a Michell 
structure within a rectangular domain. 


CONSIDERATION OF TWO SPECIFIC ORTHOGONAL COORDINATE SYSTEMS 
FOR THE MEMBER LAYOUT OF MICHELL OPTIMUM STRUCTURES 


Two specific orthogonal coordinate systems - the dipolar and the 
spherical - are selected for the member layout of Michell optimum structures 
for each system, and the strain-displacement equations (20) are solved for 
the displacements, with the strains ecxaj epp, and e a(3 corresponding to the 
Michell strain criterion. The system of curves is enclosed within some 
regular boundaries and a self -equilibrating system of forces is made to act 
along the boundary. A Michell optimum structure then obtains for the selected 
force system. The volume of this structure can now be obtained from equa- 
tion (l) since the displacements of the points of application of the forces 
are determined. Examples of some force systems are presented as applications 
of the coordinate systems for member layouts of the Michell optimum 
frameworks . 


The Dipolar Coordinate System 

The dipolar coordinate system, figure 6, is defined by the following 
transformation: 


x 


a sin P 

cosh cl - cos 3 9 


y - 


a sinh a 
cosh a - cos 3 


m 


where (a,p) represent curvilinear coordinates in the plane of the rectangular 
coordinate system and a is a constant parameter. The corresponding unit 
arc lengths A and B are determined from equations (13) and are obtained as 

A = B = — - (45) 

cosh a - cos P 


For the above values of A and B, the equation of compatibility of strains, 
equation (21), can only be satisfied if (eaa “ epp) is zero. It follows that 
this system only admits of Michell strains of the same sign. The correspond- 
ing Michell criterion for strains is taken as 

e coa ” e pp “ ±e > e ap = 0 (^6) 


Since this system can only admit Michell strains of the same sign, the 
resulting layout of the members will give an all -tension or an all- 
compression structure. The coordinate curves, along which the members will 
lie, dictate the natural boundaries and the force system that can be equi- 
librated. Inversely then, a Michell structure is obtained for that force 
system. The displacements of the system (obtained by solving the strain- 
displacement equations for the strain system of eq. (46)) give the displace- 
ments of the points of application of the forces from which the weight of the 
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structure is determined. Hence, the strain-displacement equations for the 
system will now be solved for the strain pattern of equations ( 4 - 6 ). 

The expressions for strains in this coordinate system are obtained by 
substituting the expressions for A and B from equation (4-5) into equa- 
tions (20). Then, by employing the Michell criterion defined by equa- 
tions (46), the differential equations for the displacements u and v along 
the a and 3 curves, respectively, are obtained as 
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cosh a - cos p du sin P 

a da a 


cosh a - cos P dv sinh a 

i 3p o~ 


v — ±e 


u - ±e > (47) 


cosh a - cos P dv sinh a cosh a - cos P du sin P 

+ v — + ■ ■ + u - 0 

a da a a ^p a 


When ^ = cosh a - cos P, the above system of equations is rewritten as 

t — v = ±ae 

da dP 


(48a) 


, dv d\|/ 
i|/ — _ — u - ±ae 
dp da 


k ( * T) + ^ ( * u) ' 0 


(48b) 

(48c) 


From equations (48a) and (48b) , 


7 - (tu) = ~ (tv) 

3a 30 


(49) 


Now equations (48c) and (49) give 


T~2 + = 0 atld 

3a 2 30 2 


& (*v) + ^ (»V) - o 


(50) 


Multiplying equation (48a) by ty, differentiating with respect to a, elim- 
inating (3/3a)(\|rv) by means of equation (48c), and employing the first of 
equations (50) yields a second-order differential equation that contains 
derivatives of u with respect to 0 only: 


t 




(5D 


Similarly, from equations (48b), (48c), and the second of equations (50), 
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(52) 


* 


Bf_ 

Ba 2 


(tv) 





Bt 

= ?ae Bp 


When the general solutions to equations (51) and (52) are obtained, the 
resulting expressions for the displacements u and v must be of the form 


u(a,p) = 


v(a,3) = 


f 1 (a)sin p + f 2 (a)(cos p - sech a) + ae tanh a 

i(a,p) 

f 3 (P)sinh a + f 4 (p)(cosh a - sec p) + ae tan P 

t(“,P) 


(53) 


for arbitrary functions f]_(a), f 2 (a>), f 3 (p), and f 4 (p) • Equations (48a-c) 
now determine the form of the arbitrary functions and yield finally: 


u(cc,3) = 


v(a,p) = 


(k sinh a + l cosh a) sin P + (+ae sinh a + n cosh a) cos p - n 

t(a,P) 

(n sin P - l cos P)sinh a + (+ae sin p - k cos P)cosh a + k 

t(a,P) 


(5b) 


where k, 1, and n are arbitrary constants. 

Equations (54) form the general solutions to the system of equa- 
tions (48a-c) or (47) • A general solution has been obtained for the displace- 
ments of the dipolar coordinate system for the Michell strain condition 
defined by equations (46) . If the system of curves is enclosed within suit- 
able boundaries, a system of forces can be equilibrated. Then the constants 
in solutions (54) can be determined from the physical conditions of the 
problem. This will now be demonstrated with respect to two applications of 
the dipolar system in the development of Michell optimum structures. 

Application 1 . - For the first application, a region bounded by two non- 
concentric circles a = and a = a 2 is taken (see fig. 7)* The loading on 
the resulting structure is taken to consist of uniform external and internal 
pressure p x and p 2 acting on the 3 boundaries a = a x , a = a 2 , respec- 
tively, as shown in the figure. Both sets of members are in compression. The 
general solutions for displacements for the pertinent strain system of 
e <xa = -e = epp, e ap = are given by 
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u ( a ^p) _ ae slnh a cos g + k sinh a sin p + Z cosh a sin |3 
^ \|r t 


+ n 


cosh a cos ^ - 1 




( 55 ) 


v(a,|3) = 


ae sin 3 cosh a k(l - cos (3 cosh a) Z cos j3 sinh a 
^ i \|A 


+ n 


sinh a sin g 
* 



J 


Figure 7-- Application 1 of the dipolar coordinate system: two nonconcentric circles under 

external and internal normal pressure - a candidate for entry body configuration. 


The boundary conditions are as follows: 

(i) At P(a = oo) for any 3 direction, u = 0. 

(ii) Along AB (i*e., for (3 = 0 or it and any a) due to the symmetry 
of line AB, the displacement v = 0. At (a, (3 = it), the condition v = 0 
gives (from eqs . (55)) 
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k = 


-Z sinh a 


(56) 


1 + cosh a 


At (a-, 3 = 0) , condition v = 0 gives 

l sinh a 

k = 

1 - cosh a 


(97) 


However , since k and Z are constants, it follows, from equations (56) and 
(57)> that 


k = Z - 0 (58) 

At (a = 00,3), the condition u = 0 gives 

n = -ae (59) 

From equations (55) through (59) it follows that the displacements correspond- 
ing to the structure and its loading are given by 


u(a,3) = — — (sinh a cos 3 - cosh a cos 3+1) 

cosh a - cos 3 

> 

v(a,3) = ; — 5 1 1 ■ ~ x (cosh a sin 3 - sinh a sin 3) 

v cosh a - cos 3 J 


(60) 


Once the displacements are known, the solution to the problem is complete, for 
its volume can be obtained by means of equation (l) . The volume is given by 
(for | ct c | =• I at I - a) 



where B, obtained by using equations (13), is given by B = (a/cosh a - cos 3) 
and u is given by equation (60) . 

Since the application of the dipolar coordinate system in obtaining the 
Michell structure for the loading illustrated in figure 7 has been presented, 
it would be interesting to consider some practical use of that load system. 

It is conceivable that such a configuration could be used for an entry body. 
When thus employed, the inner circle could represent the payload compartment 
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which is displaced with respect to the geometric center of the body; the 
external and internal pressures, then, would represent the aerodynamic forces 
and the inertia forces, respectively. The a members form the structural 
members to carry the forces due to the aerodynamic loading. Also, at impact 
the a members could buckle and thus act as energy absorbers. 

Application 2 (see fig. 8) .- Herein the Michell structure is bounded by 
a curves P = (B 0 an( ^ P - Jt (or the y axis) . The loading consists of a 
distributed pressure (p x ) loading along PQ (p - it) and reacted by normal 


P| applied pressure 



X 

Figure 8.- Application 2 of the dipolar coordinate system: transmission of applied pressure 

from one surface to another of greater surface area. 

pressure (p 2 ) along the surface P = P 0 . All the members, that is, the a 
and the P sets, are in compression. The Michell strain criterion is again 
e aa - -e = epp, e a p - 0 and the corresponding general expressions for the 
displacements u and v are given by equations (55) • The boundary conditions 
are as follows: 

(i) Along a = 0, for all p (i.e., along the x axis) u = 0 due to 
symmetry . 

(ii) Along P = tc (i.e., the y axis), the displacement v is zero. 
These boundary conditions give 


k = l = n = 0 


(62) 


Thus, the displacements are given by 
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(63) 


u(a,0) = — sinh a cos P ] 


* 


/ ae sin 0 cosh a 

v(a,0) = 


¥ 


The volume of the structure is now given by (for |cr c | = |a = a) 






Pi dsi J 


3 =* + P 2 ds x l 


P=3r 


2 

a 


v a=oo 


a=o 


pi i; + f p 2 i 

y 3=rt ^ e/ 3=3 0 J 


da 


( 64 ) 


where A - (a/cosh a - cos 0 ) and. v is obtained from equations (63). 

It is conceivable that this form of configuration could be used as a 
bearing structure for a load acting on a weak material or soil (i.e., one 
having a low bearing value). The shape of the outer boundary (0 - 0 O ) along 
which the material or soil forces act is obtained from the geometry of the 
dipolar system for the corresponding applied pressure distribution along 
0 — it • Along this boundary (0 =■ 0 O ) , no shearing forces are developed; con- 
sequently the possibility of slip, and hence large settlement is reduced. 

The structure then efficiently distributes this load from the bearing area 
PQ over a larger area (along the surface of the 0 = Po curve) and hence 
reduces the tendency of large deformations in the material. 


The Spherical Coordinate System for Member Layout 
The spherical coordinate system given by the transformation 

x = a sin 3 cos 7 
y = a sin 3 sin 7 \ 
z = a cos 3 


(65) 


where the curvilinear coordinates a, 0, 7 are illustrated in figure 9> is 
proposed for the member layout of a Michell optimum framework. Expressions 
for displacements are obtained for a Michell strain criterion; they represent 
the general solution of the strain-displacement equation for a Michell strain 
criterion in a three-dimensional spherical coordinate system. Any self- 
equilibrating force system can now be made to act upon a domain enclosing a 
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spherical system of curves , the proper boundary conditions on the displace- 
ments could be imposed, and the volume of the structure could then be 
det er mined - 


z 



Figure 9*~ Spherical coordinate system. 


The parameters of the coordinate system are given by 



with hg and h 3 similarly given by differentiation with respect to J3 and 7, 
respectively. For the system, represented by equations (65), the parameters 
are obtained as follows: 


hi = 1 , 



h3 


1 

a sin (3 


(67) 


The corresponding strain-displacement equations are given by (see p. 54, 
ref. 11) 
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8% d 

ea “ = hl ~ + hlh£U P Sp 


' + hlh3 ^ k 

Suq S 

e 33 = hs ~ + hsh3U 7 §7 

to 

1 + h i h 2%. ST 0 

du? d . 

e 77 = 113 “ + hlh3 ““ S 1 

(k) 

1 . . 3 /l\ 

1 + kzh3U i3 5F W 


Sex, 

i du P 




du-y 


COS P 

+ — Uo 

a sin 3 P 
hp d 


> (68) 


ii3 5 


COS 3 


l 


Sup 


e 07 * if ^ (h3U ’ ,) + c S7 (h2U P ) * - a sin e "r + a W * a sin |3 V 

. = £a a. (hl% ) + Sl 3 {h3V , ) + 

7 a hi ay ^ h3 Sa ' a sin P S7 Sa a * 


e »p * §: Jr (hsU » ) + £r Jr (hi%) ‘ 


du. 


3 _ l 


1 du** 


hi S3 


da a U P a S3 


J 


where Ua, Up, are the displacements along the a, 3, and 7 directions, 

respectively. 


For a Michell strain criterion it is necessary that 


e oa > e pp ~ e > e 77 ~® > ®ap ® 


e 3l 


(69) 


The equations (68) are now solved for the following general strain system: 

e aa > e 3p > e 77 > e ap = 6 — e p7 — e^a (76) 


It will, however, be shown that equations (68) are only solvable for certain 
relations between eoa> e 33> an( ^ <277 • These relations will then eliminate 
some of the eight cases shown in equations (69) and will be the conditions of 
compatibility of strain. For the strain system of equation (70) , the strain- 
displacement equations (68) are as follows: 


Su tx 

da e<xa 


(71a) 


26 



1 

a dp 


% _ e 

- - e pp 


( 71 b) 


du 


a sin 0 by 


+ ± 


I % + U P ■ e rr 


-cot P 

a 


1 ^ u 7 

Uy + - — '- + 


bvLr. 


P = e =0 

a d0 a sin 0 d/ ^ 


K , du 7 u 7 _ 


a sin 0 d/ da 


a 


e /a ^ 


da a P a dp 


"a0 


( 7 lc) 

( 711 ) 

( 71 e) 

( 71 f) 


Equations ( 71 a), ( 71 b), and ( 71 f ) contain only displacements u^ and Up 
and yield the general integrals 


Uq^c^P,/) = e^a + Ci(7)cos p - 0^(7) sin P ( 72 ) 

Up(a,P,7) = (epp - e aa )ap - Ci(7)sin P - Cg(7)cos P - aZ»(7) (73) 

where 01(7), C 2 (7), and Z(/) are arbitrary functions of 7 , and the primes 
indicate differentiation with respect to 7. When the derivatives of Uy 
have been eliminated between equations ( 71d) and (71e), it is seen that the 
function Cx(y) is, in fact, a constant. Equation ( 71 °) can now be integrated 
to yield 

Uy (a,P,7) = [(e 77 - e^) sin P - (e^ - epp)P cos P]a7 

+ C 2 (7) + a,Z(7)cos P + aA(a,p) (7^) 

with A(a,p) arbitrary. Now equation ( 71 e) shows that 

C 2 ( 7 ) = bi sin 7 + D 2 cos 7 + D 3 (75) 

and 

A(a,p) = ak(P) - D 3 (76) 
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where Di, D 2 , and D 3 are constants, and k is an arbitrary function of 3 . 
Finally, equation (71d) shows that 

l(y) = D 4 sin 7 + D 5 cos 7 + D 6 ( 77 ) 

where D 4 , D 5 , D 6 are arbitrary constants. Moreover, equation (71d) can be 
satisfied if and only if 

e 33 " e oa = 0 

which is the compatibility condition for the differential system (71) • With 
the arbitrary constants renamed, the general solution representing the 
displacements is written: 

u a (a, 3 , 7 ) = e nn a - A cos 3 + (B cos 7 + C sin 7 )sin 3 (79a) 



Up(a,3 , 7 ) = A sin 3 + (C cos 3 + Ea)sin 7 + (B cos 3 - Da)cos 7 ( 79 b) 


u^(a,3,7) = a, 7 (e^ - eQ^Jsin 3 + Fa sin 3 + (Da cos 3 - B)sin 7 

+ (Ea cos 3 + C)cos 7 (79c) 

From equation ( 78 ) it is seen that e rtr». and epp have to be equal in magnitude 
as well as in sign eliminating four of the eight cases of equation ( 69 ) . The 
four possible cases are as follows: 


m 


e 77 

-J-e 

+e 

•+e 

+e 

+e 

-e 

-e 

-e 

+e 

-e 

-e 

-e 


(80) 


The term eyy - e^ in equation (79c) provides for handling any of the above 
four cases. 

The formal solution of the spherical coordinate system is now complete; 
for once the displacements are known, and the volume can be obtained for a 
set of self-equilibrating external forces . 

The solution given by equations (79) for considering applications to 
some load systems shows that when eyy = e^, equations ( 79 ) can represent 
the solutions for a uniform radially loaded spherical structure. When 
eyy ~ -eoo,, U 7 - (see eq. (79°)) i s a- multivalued expression that leads to an 
interesting application. Consider a structure with a spherical boundary 
surface . Let cuts be made along two longitudinal planes , separated by a 
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small angle d/, up to the axis of the sphere. If the portion of the sphere 
enclosed between these two planes is removed and the cut surfaces are 
rejoined, a structure with initial stresses is obtained. Thus, when 
eyy = -e^, equations ( 79 ) can represent the solution for a structure with a 
spherical bounding surface with initial locked stresses in it (i.e., a 
prestressed spherical structure) and under a suitable force system. 


CONCLUDING REMARKS 


The Michell strain criterion has been successfully applied to develop 
minimum volume structures for a variety of domains to equilibrate some 
interesting force systems of significance to the aeronautics and space 
industry. 

The indeterminacy of Michell structures was discussed in the first sec- 
tion. The utility of the superposition theorem for Michell structures was 
presented, and the theorem was proved and applied to the two load systems of 
figures 1 and 2. The superposition theorem now makes it possible to solve 
a broader spectrum of force systems which by themselves do not have a unique 
solution (i.e., do not lend themselves to a unique optimum structure) but, on 
resolution into two or more systems , do lend themselves to separate solutions 
which can then be superposed. 

In the second section the governing equations, developed by Hemp, for 
the coordinate curves of the member layout of Michell structures were pre- 
sented and two solutions were obtained. These solutions give rise to the 
coordinate systems shown in figures 3 and 5 for the member layouts of Michell 
structures. These coordinate systems are then applied to offer solutions to 
a variety of force systems acting on a plane rectangular domain as well as on 
a curvilinear domain as illustrated in figure 4. 

In the third section, two well-known orthogonal coordinate systems, 
namely the dipolar and the spherical systems, were considered for member 
layout of Michell structures. General solutions were obtained for displace- 
ments to satisfy Michell strain criterion in the two coordinate systems. Two 
interesting structures arise as a result of the application of the dipolar 
system; one of them suggests interesting possibilities for an entry body con- 
figuration having a displaced center of gravity with respect to the geometric 
center, which could be helpful in orienting the body during descent. The 
other structure enables a concentrated load to be distributed over a wide 
bearing area, thereby reducing the possibility of large deformations of the 
bearing material. 

For the spherical coordinate system, the general solution of principal 
strain-displacement equations has been derived. Applications to physical 
problems could further the utility of this solution. It is shown that the 
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spherical coordinate system admits solution to the Michell strain criterion 
for the four cases listed in equation ( 80 ) . Some interesting applications of 
the solution are discussed. 


Ames Research Center 
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